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A NEW IMAGE ENHANCEMENT ALGORITHM WITH APPLICATIONS 
TO FORESTRY STAND MAPPING 

■By Edwin. P. P. Kan and Jinn~Kai Lo 
Aerospace, Systems Division,/ Lockheed Electronics Company/ Inc. 

1 . 0 INTRODUCTION 


•Satellite and aircraft multispectral scanner (MSS) data 
have been classified and analyzed for land use planning/ 
forestry mapping/ resource inventories/ and agricultural 
applications (refs. 1 through 3) . Computer classification 
maps from these studies always appear to need refinement/ 
in the sense that "holes" exist on the maps. For example, 
a geographic area supposedly containing one feature, such 
as a forest, often is. not classified entirely as a forest 
on the map. These spotty, discontinuous maps normally are 
not satisfactory to- the user, who desires smooth, clean 
maps. This report describes a new image enhancement algo- 
rithm to "clean up" computer classification maps with appli- 
- cations tp forestry stand mapping. . • . 

Forestry classification maps from the processing of 
aircraft MSS data particularly need cleaning. For example, 

with data, at a spatial resolution of less than (10 me.ters)^, 

. . 2 
i.e..., with one picture element (pixel) covering (10 meters) 

on the ground, more than one type of pixel wpuld be observed 

from aircraft over ,.a "homogeneous" timber stand. Some pixels 

■in a -Stand would represent crowns of the predominant tree 

type. Other pixe.ls would represent understories , clearances, 

and .shadowsr among ^.the... trees. A computer classification of 

such MSS. data would, .produce heterogeneous maps in which 

pixels within a stand would be classi-fied into more than 
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one statistical class. A widely accepted classification 
rule is the training— field, maxiiQum likelihood classifier 
(ref. 4) . 

The problem of holes, in classification maps is further 
compounded by the U,S. Forest Service definition of a stand 
type. A stand type refers to an area larger than 4 hectares 
in which more than 51 percent of the trees in the canopy belongs 
to one tree species or group of species. A stand smaller 
than 4 hectares should be absorbed into and classified the 
same as the surrounding area under the present system of stand 
type mapping. Because existing statistical classification 
techniques have not been designed to clean up these small 
areas, computer generated classification maps would contain 
undesirable holes which should not exist on final stand type 
maps . ■ • 


To obtain usable timber stand classification maps, the 
computer classification maps need to be postprocessed to 
remove the holes. Some investigators have used' spatial infor- 
mation to fill the holes by checking neighboring pixels as in 
the RECLASS computer program (ref. 5). In that program, a 
pixel’ surrounded by four neighboring pixels belonging to one 
different class is reclassified into the class of its neighbors. 
Variations of this neighbor rule have been explored. 

A natural generalization of the neighbor rule is that any 
connected set of pixels is eliminated if the set contains less 
than a prespecified mamber of pixels. The whole set -is then 
reclassified into the class of the surrounding pixels. By a 
connected set, it is meant that every pixel in the set has at 
least one neighbor in the same set; set with one pixel is also 
connected. This postprocessing concept is indeed appropriate 
for timber stand mapping. It will be explored in this* paper, 
formalized into a new image enhancement algorithm, and applied 
to a forestry application using aircraft MSS data. 



2 . 0 CONNECTIVITY 


A set of pixels is connected if every pixel is a left- 
right or top-bottom neighbor of at least one other pixel. 

The set with one pixel is also connected. The left-right 
connection is along a scan line, or along the same row of 
the two-dimensional array. The top-bottom connection is 
across a scan line, or along the same column of the two- 
dimens ional array . 

The five pixels labeled ”1” in figure 1(a) are con- 
nected according to the above definition, but the five 
pixels labeled "1” in figure 1(b) are not. The new image 
enhancement algorithm is developed, programed, and tested 
according to the left-right, top-bottom definition of con- 
nectivity. However, the algorithm can be generalized 
easily for diagonal connectivity in which pixels lie 

I 

diagonally adjacent to one another; figure 1(b) is an 
extreme example of diagonal connectivity. 

Formal mathematical definitions of connectivity and 
related concepts use terms such as "four-path" and "four- 
connected set." These concepts were studied by researchers 
in the pattern recognition of digital pictures, particularly 
handwritten character recognition. References 6 through 11 
cover the locating of edges in digitized pictures, thinning, 
and seeking skeletons in pictures. The new image enhancement 
algorithm is a variation of the previous research concepts 
and applies to forestry timber stand mapping. 
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3.0 THE NEW IMAGE ENHANCEMENT ALGORITHM 
3.1 CLASSIFICATION IMAGES 

The new algorithm is designed to refine classification 
images obtained by processing MSS data through techniques 
such as the training-field, maximum likelihood classifica- 
tion (ref. 4) . Classification images are two-dimensional 
arrays in which each pixel has a label representing a 
class, such as a timber type. 

To simplify the notations without loss of- generality, 
a classification image is assumed to be binary, in which 
each pixel has the labeled value 0 or 1. This two-class 
classification image is a general representation, because 
any map can be examined one feature at a time. The pixels 
labeled 1 represent the features of interest; and the 
j pixels labeled 0 are of no interest. This report presents 
the new algorithm in terms of images of 0*s and I's. 

The algorithm modifies the input binary image of 0 ’ s 
and I's. The first steps of the modification are to find, 
all the connected sets in the image and to count the nimiber 
of pixels in the sets. The labels of the sets smaller than 
a prespecified size n^ are then changed so that 1 becomes 0 
and 0 becomes 1. In this change, sets smaller than n^ 
are absorbed by the surrounding sets . 

The ways of changing labels are easy to modify. For 
example, the programed algorithm can be modified to permit 
small connected sets of I's to be absorbed and changed to O's 
and to leave small sets of O's unchanged. In section 4.2, 
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REPRODUCBILITY OF THE 
ORieiNAL PAOE IS POOR 


the tertiary images to be refined have three labels, -1, 0, 
and 1. The program was developed to modify only the two 
classes 0 and 1 and to leave the class -1 unchanged. 

3.2 SKELETON FLOW OP THE ALGORITHM 

Because of the finiteness of computer core memories, the 
image must be analyzed block by block. To comply with the 
normal line-by-line format of recording an image, a fixed 
number of lines in the image is analyzed at one time. To 
accommodate the worst case of a connected set one pixel 
wide down a column of the image, n^ lines must be studied 
at one time; n^ is the minimum size of a set not to be 
modified. After these n^ lines are analyzed, the start ’ 
line of the n^ lines is purged. The line after the stop 
line of the n^ lines is then input to the computer memory 
to make a new set of n^ lines. 

For a data image array of the size M by N in which 
the elements are 0 and 1, three integer arrays are programed 
in the computer. 

1. Data array A(I , J) : I = 1, * • ♦ ,n^ ; J = 1, • • • ,N 

2. Tag array S (I,K,L) : I = 1, • • • ,n^ ; K = 1, * • • ,N 

L = 1,2,3 

3. Count array C{P,Q); P =!,*•*, N ; Q = l,***,2nQ+3 

The use of the arrays and the skeleton of the algorithm 
are described as follows. 

1. Initialization: Input n^ lines of data and store 

them in the data array A(I,J) . 
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2. Tagging: Follow the lines of A(I,J) sequentially to 

search for along-the-line connected sets. Tag the sets 
using integers starting from 1 and store them in the 
tag array S(I,K,L). 

3. Clustering: Check the across -line connectivity of 

tagged sets. Modify the tags so that each connected 
set has- a unique tag. 

4. Counting; Count the pixels in the clustered, tagged 
sets. Store -the totals in -the count array C(P,Q) . 

5. Label Modification: Change the 0 labels to 1 and 1 to 

0 in all clustered, tagged sets smaller than the minimal 
size n^^ . 

6. Reiteration; Input a new line into A(I,J) and return to 
step 1. Stop if all lines in the data image are exhausted. 

The appendix has a detailed description of steps 1 through 6. 

3.3 TAG AND COUNT ARRAYS 

In -fehe tagging step 2 of section 3.2, the program goes 
through the data array A(I,J) . For example, in line I* 

A. 

the program finds and establishes the Kth connected set, 
stores the start pixel in S(I,K,1), stores the stop pixel 
in S(I,K,2), and tags it in S(I,K,3). Evidently, the K + 1 
connected set on the same 1th line has 


S(I, K + 1, 1) = S(I,K,2) + 1 


/s. 

*I denotes an integer variable; I • denotes a partic- 
ular integer value of I . The same applies for K, L, P, 
and Q . 
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and 


S(I, K + 1, 3) = S(I,K,3) + 1 

The count array provides N counters , each with 
2n^ + 3 entries. The N counters are sufficient for 
analyzing any image, although the image might have more than 
N connected sets. Actually the N counters are used 
cyclically. The number of times that a particular "Sth 
counter is used is stored in C(P,3). 

In the counting step 4 of section 3.2, the program,- 
goes through the tag array S(I,K,L). For example, in 

/N ^ 

line I at the Yith clustered, tagged set, the program has 
two alternatives. 

/s. /s 

1. The Kth set is connected to the above line I - 1 . 

/V , /\ 

2. The Kth set is not connected to the above line 1^1, 
and is a new clustered, tagged set. 

In alternative 1, the program looks up the Yth counter, 
where P = mod[S (I,K, 3) ,N] . The mod(M,N) denotes the 

/N. 

remainder of M divided by N . The C(P,1) contains a 

^ A. 

Q value, where 4£Q£2n+ 3 , such that the program 

^ A /V 

stores the start pixel of the Kth set in C{P,Q) and the stop 

A A 

pixel of the set in C(P, Q + 1). The program increases the 

A 

current Q by 2. In alternative 2, a new counter is used 

A 

for storing the start and stop pixels of the Kth set. The 
entries C(P,4) to C(P, 2n^ + 3) of any Vth counter store 
the start and stop pixels of the clustered, tagged sets 
from the first line until the last line in which these 
sets appear. The C,(P,2) stores the current total number 
of pixels in the clustered, tagged sets. 
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4..-0' TIMBER STAND MAPPING 

A goal of the Forestry Applications Exploratory 
Studies Project (ref. 12) is to map timber stands by 
computer-aided analyses of MSS data collected by aircraft 
and satellites. The MSS data collected over the Sam Houston 
•National Forest in east Texas were analyzed with a computer 
algorithm of the LARSYS-type (ref. 4) . Computer classifi- 
cation maps were generated and postprocessed by the new 
image enhancement algorithm. 

4.1 DATA SETS 

The three MSS data sets used are IX, 2X, and 3X. The 
IX data were collected during Mission 230 of the C-130 air- 
craft of the National Aeronautics and Space Administration 

on March 21, 1973, at a 3-kilometer altitude. The data have 

2 

a spatial resolution of (8 meters) ; i.e. , a pixel covers an 
8- by 8-meter area on the ground. In the Forestry Applica- 
tions Exploratory Studies Project, the IX data set is called 
edit Si (ref. 13) and covers about 250 scan lines with' 

700 pixels per line. 

The 2X and 3X data are electronically reduced from IX 

by pixel- averaging (ref. 14) . The 2X and 3X data cover the 

same geographic area of IX but have spatial resolutions of 
2 2 

(16 meters-) and (24 meters) , respectively. The numbers 
of scan' lines and pixels per line are reduced correspondingly 
by two and three times. The 2X and 3X data simulate the data 
that would be acquired at two and three times the altitude of 
IX. 
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Figure 2 is a three-channel rendition of the IX, 2X, 
and 3X MSS data. Because of the pixel-averaging process 
producing 2X and 3X data from IX, the 2X data have a smoother 
texture than IX, and the 3X data have a smoother texture than 
the 2X. The physical ground features on the data sets consist 
of a regeneration stand, pine stands, mixed hardwood stands, 
cutover land, and open areas of roads, pipelines, and ponds. 

Using the same parameters in a clustering procedure, 
the -iX, 2X, and 3X data were analyzed, resulting in the 
cluster classification maps in figures 3(a)_, 4(a), and 5(a). 

A clustering procedure groups the pixels with similar data 
values into clusters, or homogeneous spectral signature 
classes. A cluster classification map is the image of the 
clustered data in which each pixel is labeled with the class 
to which it belongs. 

Figures 3(a), 4(a), and 5(a) show three-class decom-. 
positions of the data. The lightest shade is pine stands; 
the medixam shade, hardwood stands, cutover land, and regen- 
eration areas; and the darkest shade, open areas. These 
figures are by no means -ultimate timber stand classification 
maps, but they suffice here to illustrate the image enhance- 
ment power of the new algorithm. 

4.2 APPLYING THE ALGORITHM 

Figures 3(a), 4(a), and 5(a) are tertiary images with 
the labels -1, 0, and 1. The algorithm changes 0 to 1 and 
1 to 0 when necessary, leaving -1 unchanged. The -1 label 
corresponds to the. open area represented by the darkest shade 
in the maps. 
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Figure 2.— 3-channel color renditions of MSS/24 data. 
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gure 5.— Forest classification maps before and after 
'ostprocessing Sam Houston National Forest 3X data. 
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The algorithm was applied to the cluster classification 
map for IX data for values of 5^ 10, and 40. Connected 

sets smaller than n^ were absorbed by the surrounding sets . 
The enhanced maps are in figures 3(b), 3(c), and 3(d). 
Similarly, figures 4(b), 4(c), 4(d), 5(b), 5(c), and 5(d) 
show the enhanced maps of 2X and 3X data for n^' values of 
5, 10, and 40. 


4.3 RESULTS 

Figures 3 through 5 show how the new image enhancement 

algorithm can remove the holes in classification images while 

preserving significant boundary information. For the IX data 
2 

at a (8 meter) spatial resolution, postprocessing with 
n^ = 40 produced the best result. The same n^ value 
produced the same result for 2X data. For the 3X data, 
postprocessing with n^ = 10 provided the best result 
because n^ = 40 produced a map in which some significant 
features were altered. These results show the possible 
compromise between data spatial resolution, classification 
accuracies, and optimal postprocessing with the best 
n^ value. Indeed, another implication from the results 
was the possible advantage to timber stand classification 
by using lower resolution data (i.e., data having large 
on-ground pixels) ; such a conclusion had been drawn from 
previous works (ref. 14) . 

The new image enhancement algorithm removed from the 
original IX data the variously shaped, connected holes 
smaller than the n^ values . Linear features were retained 
except in figure 5(d) , in which n^ = 40, was too large for 
data at that spatial resolution. The curvilinear shape of 
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the boundaries was preserved, not smeared. These set- 
removal, shape-preservation properties are needed in timber 
stand mapping. These properties are weak or absent in 
other postprocessing techniques, such as neighbor-checking. 
Section 5.0 demonstrates via examples the superiority of the 
new image enhancement technique, in regards to such properties 
which are required in timber stand mapping. 




?AQfS ms 







5.0 THE NSW IMAGE ENHAHCEMEMT ALGORITHM AMD 
NEIGHBOR-CHECKING PROCEDURES 


5.1 THE NEIGHBOR-CHECKING ALGORITHM 


A neighbor-checking algorithm, such as RECLASS, post- 
processes a classification map by modifying the pixel labels 
when a majority of the surrounding pixels have different 
labels. Sections 5.1 through 5.5 refer to the eight- 
neighbor- checking algorithm in which m^ is the average 
of all eight neighboring labels surrounding x (fig. 6) and 
T is the input threshold in which 0.5 £ T £ 1 . The 
algorithm changes x to 1 if and only if T _< m 5 1 • The 
algorithm changes x to 0 if and only if 0 £ < 1 - T . 

Neighbor checking always uses the label values of the 
map before relabeling, which retains the old values until 
the checking is complete. 

For T = 0.5 , the center pixels x^, x^, x^, x^, and 
Xg , respectively, in figures 6(a), 6(b), 6(c), 6(d), and 
6(e) will be changed to I’s by this neighbor-checking algo- 
rithm. Similarly, x^, x_, x*, and x* , respectively, in 
figures 6-(f ) , 6(g), 6(h), and 6 (i) will be changed to O’s. 
The modification is irrespective of the original values of 
the center pixels. 


5.2 SET REMOVAL 


The neW' image enhancement algorithm can remove unwanted 
isolated sets, i.e. , isolated areas (fig. 7) , more easily 
than the neighbor-checking algorithm. Isolated ^^^eas in 
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Figure 7.— An isolated set of I’s removable by the new 

algorithm with n^ ^ 13 but unremovable by the neighbor- 
checking algorithm with 0 . 5 ^ T ^ 1 . 
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timber classification maps must be removed when the area size 
is smaller than n^ . 

The isolated set of I's in figure 7 is removable by the 
new algorithm if > 13 . However, using 0.5 £ T £ 1 in 

the neighbor-checking algorithm, the eight neighbors of each 1 
have the mean value m > 0.5 . The neighbor-checking algo- 

X — 

rithm can remove none of the I's in the figure. 

•An "unremovable set" is a connected set which cannot be 
eliminated by label modification by applying a given algo- 
rithm a finite number of times. The isolated set in figure 7 
is. removable by the new algorithm but unremovable by' the 
neighbor-checking algorithm. This fact holds for any con- 
nected set containing the isolated set in figure 7 . 

In the smaller isolated set in figure 8(a) the I's are 
removable by the new algorithm with n^ ^ 6 . The first use 
of- the neighbor-checking algorithm with T = 0.5 reduces 
the set size but fails to remove the set. See figure 8(b). ' 
The second use of the neighbor-checking algorithm removes 
the set. 

Therefore, the isolated set in figure 8(a) is removable 
by. the new algorithm with n^ ^ 6 . The set is removable by 
the neighbor-checking algorithm on the second try. 

In conclusion, isolated data sets as in figures 7 and 8 
must be removed for timber mapping. The new image enhance- 
ment algorithm is more powerful in removing isolated sets 
than the neighbor-checking procedure. 
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(a) Original isolated set 
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(b) Modified image after one use of the algorithm 


Figure 8.— An isolated set of I's removable by two uses of 
the neighbor-checking algorithm with T = 0.5 ; the 
resulting modified image after two uses contains all 0 ' s 
and is not shown here. This set of I's is removable by 
the new image enhancement algorithm with ^ 6 . 
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5 . 3 LINEAIl FEATURES 

The new image enhancement algorithm preserves linear 
features such as roads and pipelines; but the neighbor- 
checking algorithm eliminates them in most cases . Linear fea- 
tures often form administrative boundaries in timber stand 
maps. The long, narrow appearance of these features is not 
just useful but sometimes mandatory, and a good post- 
processing image enhancement technique should preserve them. 
The data sets in this report all have these features . 

For example, figure 9 shows two linear features rep- 
resented by I's. The new image enhancement algorithm with 
n^ ^ 8 preserves the two linear features. However, the 
neighbor-checking technique with 0 • 5 _< T _< 1 eliminates 
the linear features completely. 


5.4 CURVILINEAR BOUNDARIES 


The new image enhancement algorithm preserves curvi- 
linear boundaries, but the neighbor -checking algorithm does 
not. Curvilinear features, as in figures 10(a) and 11(a), 
often appear in timber stand maps. Preservation of these 
features during postprocessing is essential. 

With the appropriate size of n^ , the new algorithm 
preserves the curved boundary between the 0 ' s and I's on 
figures 10(a) and 11(a). The size of n^ depends on whether 
figures 10(a) and 11(a) are parts of larger images or stand- 
alone images. 
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Figure 9.— Linear features shown by two sets of I's. The 
new algorithm with ^ 8 preserves the features, but 

the neighbor-checking algorithm with 0.5 < T < 1 removes 
them . ~ 
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(a) Original image 
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(b) First use 



(c) Second use 
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(d) Third use 


(e) Fourth use 


Figure 10.— The change in the curvilinear boundary between 
the 0 ' s and 1 ' s after uses of the neighbor-checking 
algorithm with T > 0 . 5 . 


































































(a) Original image 
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(c) Second use 


Figure 11,— The change in the curvilinear boundary between 
the O's and I's after uses of the neighbor- checking 
algorithm with T ^ 0.5 . 
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However, the first use of the neighbor- checking algo- 
rithm with T > 0,5 on figure 10(a) changes the boundary 
to that of figure 10(b). Three . successive applications 
reduce the image to that in figures 10(c), 10(d), and 10 (e) . 
The protruded portion of I’s eventually becomes smeared and 
removed, as in figure 10(e). Similarly, two uses of the 
neighbor- checking algorithm with T ^ 0.5 on figure 11(a) ' 

remove the protruded portion of I's, as shown in figures 11(b) 
and 11(c) . 


5.5 STRUCTURAL PATTERN RECOGNITION 

In structural pattern recognition, the new image 
enhancement algorithm treats classification maps on a higher 
level than the localized neighbor-checking procedure. Struc- 
tural pattern recognition is the algebraic or syntactic 
approach (ref. 15) . 

The new image enhancement algorithm decides on modifi- 
cation or nonmodification until a connected set is detected. 
However, the neighbor-checking procedure localizes its 
decision based on the labels of pixels surrounding the 
center pixel. Both the high-level and low-level treatments 
of the structural approach have advantages and disadvantages, 
the discussion of which lies beyond the scope of this report. 
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6.0 CONCLUSIONS 


The report shows the theory and applications of a new 
image enhancement algorithm which postprocesses and refines 
computer classification maps of multispectral data. The 
development of the procedure resulted from the need to 
remove holes in the data to obtain more useful timber stand 
classification maps by computer processing. The new algo- 
rithm was applied to three data sets collected over the 
Sam Houston National Forest at three different spatial data 
resolutions . 

The results show that the new algorithm produced cleaner 
classification maps' in which holes of small predesignated 
sizes were eliminated and significant boundary information 
was preserved. These cleaner postprocessed maps better 
resemble true-life timber stand maps and are thus more usable 
products than the pre-postprocessed ones. Compared to an 
accepted neighbor-checking postprocessing technique, the new 
algorithm is indeed more appropriate for timber stand mapping 

More work remains to be done, as in the relationship 
between spatial data resolution and classification accuracies 
optimal usage of the procedure, and optimal parameters used 
in the algorithm. Variations of the new image enhancement 
should be applied to more, varied data sets. The computer 
time required to perform the algorithm should be compared to 
the time required by other algorithms. The overall effect 
of treating multiclass images as binary images of 0 ' s and 
1 ' s should be examined as to the modification of connected 
sets during step-by-step applications of the algorithm to 
individual features. 
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7.0 • SUMMARY 


The theory and applications are presented of a new image 
enhancement algorithm which refines computer classification 
maps of multispectral data. The refinement eliminates 
connected sets smaller than a prespecified size and merges 
them to the surrounding area. Such a postprocessing is 
motivated by and accomplishes the requirement of a forestry 
application of timber stand mapping, where small areas are 
conventionally required to be absorbed into surrounding 
large areas to form "homogeneous" stands. The computer 
algorithm was developed, programmed, and applied successfully 
to three rauitispectfal data sets collected over Sam Houston 
National Forest at three different spatial data resolutions. 

A comparison against an accepted neighbor-checking post- 
processing technique reinforced the conclusion on the appro- 
priateness of the new procedure for timber stand mapping. 


This work was motivated by the necessity of cleaning 
"holes” in computer generated classification maps of forestry 
multispectral data. An area that is conventionally called a 
stand of forest often fails to be entirely classified by 
computers as forest due to the norihomogeneity within the 
stand, or due to openings between tree crowns, or even due 
to poor classification designs. The problem of holes is 
further compounded by the U.S. Forest Service definition of 
a stand type. A stand type refers to an area larger than 
4 hectares in which more than 51 percent of the trees in the 
canopy belongs to one tree species or group of species. A 
stand smaller than 4 hectares should be absorbed into and 
classified the same as the surrounding under the present 
system of stand type mapping. Because existing statistical 
classification techniques have not been designed to clean 
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up these small areas, computer generated classification maps 
would contain undesirable holes which should not exist on 
final stand type maps. 

A new computer algorithm was developed and programmed 
to solve the above problem with the idea of searching for 
connected sets, determining their sizes, and eliminating 
those small sets. The minimum size, , of a connected set 
not to be eliminated can be specified by the user- A set of 
pixels is connected if every pixel is a left-right or top- 
bottom neighbor of at least one other pixel -in the set; the 
set of one pixel is also connected. Diagonal connectivity 
can be similarly examined but not programmed for the present 
study. This idea of connectivity was studied by past 
researchers in the pattern recognition of handwritten 
characters. However, the present variation of the connec- 
tivity concept and its application to timber stand type 
mapping are novel. 

Classification images considered in the study were binary 
images, in which each pixel has the label value O or 1. This' 
two-class classification image is a general representation, 
because any map can be examined one feature at a time. 
Modification of the computer program to multi-class images 
is relatively simple. A typical classification image is one 
derived from processing multispectral scanner data through 
techniques such as the training-field, maximvim likelihood 
classifier. 

The new image enhancement algorithm was applied to post- 
process classification maps of the Sam Houston National Forest 
of Texas. The three data sets used were IX, 2X, and 3X, 
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* , 2 2 
having spatial resolution of (8 meters) , (16 meters) , and 

2 

(24 meters) , respectively. 2X and 3X were electronically 
reduced from IX by pixel-averaging. 

Results of applying the algorithm showed that cleaner 
classification maps were produced, in which holes of small 
predesignated sizes were eliminated and significant boundary 
information was preserved. These cleaner postprocessed maps 
better resemble true-life timber stand maps and are thus more 
usable products than the pre-postprocessed ones. 

Comparison was made via examples to demonstrate the 
appropriateness of the new algorithm against another widely 
accepted postprocessing algorithm, the neighbor-checking 
algorithm. Superiority of the new algorithm was demonstrated 
in regards to the properties of set-removal, shape -preserva- 
tion of linear features and of curvilinear boundaries. These 
properties are required in timber stand mapping. It was 
shown that the new algorithm can remove unwanted isolated 
sets more easily. Important linear features in timber type 
maps such as those due to roads and pipelines, are preservable 
by the new algorithm while usually eliminated by the neighbor- 
checking algorithm. Curvilinear boundaries are preserved 
fay the new algorithm instead of being smeared and distorted 
by the other procedure. 

In conclusion, the new image enhancement procedure is 
useful in post-processing computer generated classification 
maps and producing more usable products that better resemble 
true-life timber stand maps. More work remains to be done 
such as the relationship between spatial data resolution and 
classification accuracies, and optimal use of the algorithm. 
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APPENDIX 

THE SIX STEPS OF THE NEW ALGORITHM 


The six steps of the new image enhancement algorithm 
are initialization, tagging, clustering, counting, label 
modification, and reiteration. The notations are 

1. Data array A(I,J) : I = l,»»»,nQ ; J =!,••*, N . 

2 . Tag array S (I ,K,L) : I = 1, • • • ,n^ ; K = 1, • • • ,N ; 

L - 1,2,3 . 

3. Count array C(P,Q): P =!,•••, N ; Q = l,***,2n^+3 . 

4. Input data image D(I,J): I = 1,***,M ; J= 1,***,N . 

The integer arrays take values of 0 and 1. The symbol 
as in X := Y , denotes the placement of the value 
of Y into X . 


A.l INITIALIZATION 

The initialization steps follow. 

1. Initialize A (I, J) , S(I,K,L), andC(P,Q). 

2. A(I,J) ;= D(I,J) . J = !,•••, ,N ; I = l,***,n^ , 

3. LL := 1 ; tag := 0 . 

4. Do the steps in section A. 2 for I = l,*»*,nQ . 

A. 2 TAGGING 

Along-the-line connected sets are found and tagged. 

1. K := 1 . 

2. Along the Tth line of A{I,J) , find the Kt?i connected 
set. Then tag := tag + 1 . 
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3. S(I,K,1) := the J coordinate along the lih line of 
the first element of the Kth connected set in A(I,J) . 

4. S(I,K,2) the J coordinate along the 1th line of 
the last element of the Kth connected set in A(I,J) . 

5. S(I,K,3) := tag ; K := K + 1 . 

6. Repeat steps 2 through 5 until all connected sets on 
the 1th line are tagged. 

A. 3 CLUSTERING 

Across-line connected sets are identified and tagged 
uniquely. 

1 » I L I — LL t 1 . 

2. PP := 1 ; R1 = modCiL - 2, n^) + 1 ; R2 = mod(IL “ 1, n^) + 1 

3. Test these conditions: 

a. S{R2,PP,1) < S(R1,K,2) . 

b. S{R2,PP,2) > S(R1,K,1) . 

c. a[R1,S (R1,K,1)] = a[r2,S(R2,PP,D] . 

4. If the conditions in step 3 hold/ 

S(I,K,3) ;= min[s (Rl/K/3) /S (R2/PP,3)] I=Rl/***/R2 

5. Check each tag of all connected sets along the Rlst and 
R2nd lines. If the tag equals max[s (R1,K,3) ,S(R2/PP,3)] / 
change it to min[s (R1,K, 3) ,S (R2 /PP ,3)] . 

6. Repeat steps 3 and 4 for K = 1,2/ •••/N . 

7. PP := PP +'l . Repeat steps 3 through 6 until PP — N . 

8. If IL = n^ , go to step 10. 

9. IL := IL + 1 . Go to step 2. 
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10. IL := IL‘- 1 . Repeat steps 2 through 7 and 10 until 
IL = 2 . 


A. 4 COUNTING 

The pixels of clustered/ tagged sets are counted and 
stored in C(P,Q). 

1. I := mod(LL - 1, n^) + 1 , 

2. Compute and for all K by the equation 

S(I,K/3)/N = Qj^-N + \ . 

3. If C(R^,1) = 0 , C(R^/2) . Put S(I,K/1) and 

Jx Jv K 

S{I,K,2) in the two next available words in the R th 

Jx 

counter. The address of this next available word is 
stored in C(R^/3) . 

K. 

4. If 0 < C(R„,1) < n^ and if = C(R^,2) , put 

S(I/K,1) and S(I,K/2) in the two next available words 
in the R th counter. 

5. I ;= mod(I + 1, n^) . 

6. Repeat steps 2 through 6 the number of times shown by 


A. 5 LABEL MODIFICATION 


The tags of clustered, tagged sets smaller than n^ 
are modified. 

1. I mod(LL - 1, n^) + 1 . 

2. In checking through the counters P =!,••♦, N , if 
0 < C(P,1) <11^ , change the values in A(I,J) for 

the J locations C(P,4) and C(P,5), between C(P,6) and 



C(P,7), and so forth. The change is A(I,J) ;= 0 if 
A(I,J) = 1 ; A(I,J) := 1 if A(I,J) ^ 0 

3. Output the Ith line of A(I,J) into the final post- 
processed image. 


A. 6 REITERATION • 

‘Input a fresh line of data from D(I,J) into A(I,J) 
until all lines are exhausted. 

1. LL ;= LL +1 if all lines are not exhausted. Stop if 
they are exhausted. 

2. R J— LL t n^ . 

3. I := mod(R - 1, n^) + 1 . 

4. A(I,J) := D(R,J) . J = 1,***,N . 

Go to the steps in section A. 2. 


5. 








